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■ Abstract: We give the distribution of M n , the maximum of a sequence of n observations from a 
moving average of order 1. Solutions are first given in terms of repeated integrals and then for the 

| case where the underlying independent random variables have an absolutely continuous density. 

I/"") ■ When the correlation is positive , 

o 

O 



P(M n < x) = ^p jx v% « B x u& 

oo ; 3=1 
o 

where {y 3X } are the eigenvalues (singular values) of a Fredholm kernel and v\ x is the eigenvalue 
of maximum magnitude. A similar result is given when the correlation is negative. The result is 
analogous to large deviations expansions for estimates, since the maximum need not be standardized 
to have a limit. 

For the continuous case the integral equations for the left and right eigenfunctions are converted 
to first order linear differential equations. The eigenvalues satisfy an equation of the form 

oo 
i=l 

for certain known weights {wi} and singular values {9i} of a given matrix. This can be solved by 
truncating the sum to an increasing number of terms. 



1 Introduction and Summary 

Little is available in the literature on the behaviour of extremes of correlated sequences apart from 
some special cases involving Gaussian processes and some weak convergence results. For example 

1 Work begun while visiting the Statistics Dept, UNC, Chapel Hill. 
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Leadbetter et al. (1983) p59 give a convergence in distribution for the scaled maximum of a 
stationary sequence, and Resnick (1987) p239 gives a similar result for moving averages. 

This paper gives a powerful new method for giving the exact distribution of extremes of n 
correlated observations as weighted sums of nth powers of associated eigenvalues. The method is 
illustrated here for a moving average of order 1. 

Let {e{\ be independent and identically distributed random variables from some distribution F 
on R. Consider the moving average of order 1, 

Xi = ej + pej_i 

where p ^ 0. In Section 2 we give expressions for the distribution of the maximum 

M n = max X: 
i=i 

in terms of repeated integrals. This is obtained via the recurrence relationship 

G n (y) = I(p < 0)G n _ 1 (oo)F(y) + £G n _i(y) (1.1) 
where G n (y) = P(M n < x, e n < y), (1.2) 

1(A) = 1 or for A true or false, and K, is an integral operator depending on x. (Dependence on 
x is suppressed.) For this to work at n = 1 we define 

M = -oo so that G (y) = F(y). 

In Section 3 we consider the case when F is absolutely continuous with density f(x) with respect to 
Lebesque measure. In this case we show that corresponding to /C is a Fredholm kernel K(y, z). We 
give a solution in terms of its eigenvalues and eigenfunctions. This leads easily to the asymptotic 
results stated in the abstract. 

Our expansions for P(M n < x) for fixed x are large deviation results. If x is replaced by x n 
such that P(M n < x n ) tends to the GEV distribution, then the expansion still holds, but not the 
asymptotic expansion in terms of a single eigenvalue, since this may approach 1 as n — > oo. 

Set J r = f r(y)dy. 
2 Solutions using repeated integrals. 



For n > 1, G n of (jl.2p satisfies 

G n (y) = P{M n -i < x, e n + /?e n _i < x, e n < y) 

= P(M n _i < x, e n -i <(x- e n )/p, e n < y) if p > 



V G n „ 1 ((x-w)/p)dF(w) 
P(M n _i < x, e n _i > (x - e n )/p, e n < y) if p < 

y [G n _i(cx)) - G n _i((x - w)/p))dF(w) = G n _i(oo)F(y) - J G n _i((z - w)/p)dF(w) 

(Thanks to the referee for noting a slip in the last line.) That is, for n > 1, (jl.ip holds with 

ry 

K,r(y) = sign(p) / r((x — w) / p)dF(w). (2.1) 
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Our goal is to determine 

u n = P{M n <x) = G n (oc). 
In this section we give u n in terms of 

v n = [K n F(y)] y=00 . (2.2) 

For example 

v x = - J F(z)dF(x - pz) = -I(p < 0) + J F(x - pz)dF(z). (2.3) 

The behaviour of u n falls into two cases. 
The case p > 0. For n > 1, 

u n = v n (2.4) 



since 



G n {y)=K n F{y). 

The marginal distribution of X\ is u\ = v\ given by (|2.3[) . 
The case p < 0. By §TQ), for n > 0, 

G n +i{y) = u n F(y) +KG n (y) 

n 

= a n (y) ®u n + a n +i(y) where aj(y) = VF(y), a n (g> fr w = ^ a.jb n -j- (2.5) 

j=0 

Putting y = oo gives the recurrence equation for ii n : 

rt 

ii = 1, n n+ i = v n+ i + ^ViUn-i, n > 0. (2.6) 

i=0 

The marginal distribution of X\ is «i = 1 + v i of (|2.3[) . 

An explicit solution for u n when p < 0. 

Define the generating functions 

oo oo 

U(t) = J2^nt n , V(t) = ^v n t n , 

n=0 n=0 

Multiplying ()2.6|) by t n and summing from n = gives (?7(t) — V"(t))/t = C/(t)F(i), so that 

oo 

U(t) = (1 - W(t)) _1 F(t) where W(t) = tV(t) = ^ w n t n , w n = v n - U 

n=l 

oo 

tu(t) = (i - w^t))- 1 ^) = (i - wit))- 1 -i = Yl w ( t ) j - 

0=1 

By definition, for j = 0, 1, • • • 

oo 
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where B n j(w) is the partial ordinary Bell polynomial in w = (wi,w 2 , • • • ) tabled on p309 of Comtet 
(1974). For example 

B n o(w) = 5 n0 , B nl (w) = w n , B nl (w) = wx. 

So 

oo 

tU(t) = (1 - Wit))' 1 -1 = J2 B n {w)t n 

n=l 

where 

n 

B n (w) = y~]B nj (w) 

3=0 

is the complete ordinary Bell polynomial. For example Bq(w) = 1. Taking the coefficient of t n gives 
the explicit solution 

u n _! = B n (w), n > 1. (2.7) 

For example 

u = B~i(w) = B\\(w) = wi = v = 1, 

ui = B 2 (w) = B 2 i(w) + B 22 (w) = w 2 + w\ = «i + 1. 

Similarly from Comtet 's table we can immediately read off u n , 1 < n < 9: 

U = 1, Ul = Vl + 1 
«2 = V2 + 2U1 + 1, 

^3 = ^3 + (2u 2 +vf) + 3ui + 1, 

u 4 = w 4 + (2w 3 + 2viv 2 ) + (3v 2 + 3u?) + 4vi + 1, 

^5 = ^5 + (2t>4 + 2fit>3 + v\) + (3f 3 + 6t>if 2 + u?) + (4t>2 + 6vf) + 5ui + 1, 

= vq + (2v 5 + 2v\ V4 + 2u 2 v 3 ) + (3w 4 + 6V1V3 + 3^2 + 3vfv 2 ) + (4w 3 + Y2v\v 2 + 4«f ) 

+(5«2 + lOuf ) + 6v x + 1, 
u 7 = v 7 + (2u 6 + 2uiu 5 + 2u 2 W4 + uf) + (3v 5 + 6V1V4 + 6v 2 V3 + 3^1^ + 3viv 2 ) 

+(4t> 4 + 12viu 3 + 6u| + 12u?«2 + v\) + (5v 3 + 20uiu 2 + 10wf ) + (6u 2 + 15wf ) + 7«i + 1, 

1*8 = u 8 + (2^7 + 2v\Vq + 2v 2 1>5 + 2W3V4) + (3^6 + + 6v 2 V 4 + 3f| + 3^U4 + &ViV 2 Vz + wf) 

+(4u 5 + 12viU4 + 12^2^3 + 12u?u 3 + 12u x ^ + 4^ 2 ) + (5v A + 20uiu 3 + 10u| + 30u?u 2 + 5«f ) 
+ (Qv 3 + 30viu 2 + 20vf) + (7v 2 + 21u?) + 8«i + 1, 
ug = vg + (2v 8 + 2v\v 7 + 2v 2 vq + 2v 3 v 5 + vl) + (3^7 + 6uiu 6 + 6v 2 v 5 + 6W3W4 + 3vfv 5 + 6viv 2 v 4 
+3vivl + 3v 2 v 3 ) + (4u 6 + 12ui v 5 + 12u2f 4 + 6u| + 12vfv4 + 24viv 2 v 3 + 4uf + 4ufu 3 
+6V1U2) + {5v 5 + 20t;it;4 + 20^2^3 + 30^3 + 30v!V 2 + 20vfv 2 + vf) + (6v 4 + 30tjiw 3 + 15v 2 
+60vfv 2 + 15vf) + (7v 3 + 42v lV2 + 3hv\) + (8v 2 + 2Sv\) + 9«i + 1. 

More generally any u n can be obtained from (|2.7I) using the recurrence relation 

b n = w n (&b n , n > 1, where u>o = 0, 6 n = B n (w). (2-8) 

For example since 60 = lj this gives 

61 = tui, b 2 = w\ + u; 2 , ^3 = ^1 + 2wiu 2 + u> 3 . 

The recurrence relation (| 2 . 8 1) for the complete ordinary Bell polynomials follows by taking the 
coefficient of t n in (1 — w)- 1 — 1 = u>(l — u>) _1 where it; = W(t), and appears to be new. 
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3 The absolutely continuous case. 



Our solutions (|2.4p . (|2.6p . (|3,8p do not tell us how u n behaves for large n. Also they require repeated 
integration. Here we give solutions that overcome these problems, using Fredholm integral theory 
given in Appendix A. Write (|2.ip in the form 



Kr(y) = J K(y, z)r(z)dz where K (y, z) = pl(x < y + pz)f(x — pz). (3-1) 

Since 

«-//*(,,,)*(.,»)** - f J //(,<, + „)/(,<, + „)/(, -,»)/(,-„)** 

< P 2 f f f{x- pz)f(x - py)dydz = 1, (3.2) 



K (y, z) is said to be a Fredholm kernel w.r.t. Lebesgue measure, allowing the Fredholm theory 
of the Appendix to be applied, in particular the functional forms of the Jordan form and singular 
value decomposition. If say, < p < 1, then one can show that 

II/CHI = [FtoUM T 1 a, x T oo where *« = mi„(. - P t, ( X - ,)/,). 



Let {Xj,rj,lj : j > 1} be its eigenvalues (singular values) and associated right and left eigen- 
functions ordered so that |A 3 -| < |A J+ i|. By Appendix A these satisfy 



\jKLvj — Tj 7 XjljK^, — lj 



, J rjlk = S jk , (3.3) 



where 5jk is the Kronecker function and we write f a(y)b(y)dy = f ah. So {Vj(y), lk(y)} are biorthog- 
onal functions with respect to Lebesgue measure. Set 

By (J32D and 

oo 

where Vj are the singular values, or if the Jordan form is diagonal, the eigenvalues. (We shall use 
these terms interchangeably.) So \i>j\ < 1 and 1 + Vj > 0. 

Consider the case where the Jordan form is diagonal. Suppose that the eigenvalue Ai of smallest 
magnitude has multiplicity M (typically 1). Set 

. M 

P j =r J (^) ]Fl i ,B = Y i P j . (3.4) 

i=i 



Then by ([A8j) for n > 1, 



Y,Prf = B^{l + e n ) (3.5) 
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where e n — > exponentially as n — ► oo. (In fact by (|A.8j) 1 = i>o = X^=i /^j if this converges.) So 
for n > 1, by ([MD 

oo 

for p>0, u n = ^/3^. (3.6) 

3=1 

v\ is given by (|A,10j) with p Lebesgue measure. (When p = then (|3.6p holds with /3j = 6ji, v\ = 
F(x). So we expect that v\ — > F(x) as p I 0.) 

Now suppose that p < 0. 

By (|33D, for max™! < 1, V(t) = 1 + Y%li Pj^jt/i 1 ~ Vjt). So 

oo 

l-tV(t) = l-t-^Tp j v j t 2 /{l-v j t) = N(t)/D(t) 

3=1 

where D(t) = U°° =1 (l - v 3 t), N(t) = H°° =1 (l - Wj t) say. 

D{t) is the Fredholm determinant of K(x, y). (Now Wj takes on a different meaning than in Section 
2.) So by the partial fraction expansion, assuming that {wj} are all different, 

oo 

iV^) -1 = y - Wjt) -1 where Cj = 11^(1 - Wk/wj), 

3=1 

oo oo 

= ^2N n t n where N n = J2 c J lw j- ( 3 - 7 ) 

n=0 j=l 

Also by Fredholm's first theorem - see for example, p47 of Pogorzelski (1966), 

D(t) = l + f^D n (-t) n /n\, D n = J ••• J N^['' S s n ^d Sl ---ds n 
where -/V(^^* n ) = det(N(sj, Sk), 1 < j, k < n). Alternatively, a simple expansion gives 

D n / n \= "h---vj n = n 

l<31<-<3n 



the augmented symmetric function, in the notation of Table 10 of Stuart and Ord (1987). This 

"voo r 
^=l V i- 



table gives [l n ] in terms of the power sums (r) = Y^T=\ u j- For example [l 3 ] = 2(3) — (2)(1) + (l) 3 . 



In our case 

oo 



K r (x, x)dx = u 



r 
3 

3=1 



where by (|A.4j) 

oo 

K r (x,y) = K7' 1 K(x,y) = Y j v r j r j {x)l j {y). 

3=1 

So [l n ] has the form 

n 

t in ] = E E A(n l ---n k ){n 1 )---{n k ). 

k=l niH \-n k =n 
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However this does not give its behaviour for large n. At any rate, we have 
1 + tU{t) = (1 - Wit))' 1 = (1 - tVit))- 1 = D(t)/N(t) 

so that 

«n~i = N n <g> D n (-l) n /n\ (3.8) 

where N n is given by (|3.7p . This solution will be useful for large n if D n has an expansion of 
the form f)3.5[) . However to date we have not been able to show this directly. One can show that 
D„ = (—l) n B n (d) where B n (d) is the complete exponential Bell polynomial, d r = — (r — l)!u>(n), 
and w(n) = Y^jLi w l = ( n ) f° r w - We conjecture that if d n = J2 < jLi a j wr j w here \wj\ is strictly 
decreasing and \aA > 1, then 

B n (d) « d n « aiiy" asm oo. 
An alternative approach is to try a solution for u n of the form (|3.5|) , say 

oo 

«n = X><5? (3.9) 

i=i 

where 5j decrease in magnitude. Assuming that {Sj,Uj} are all distinct, substitution into the 
recurrence relation (|2.6p gives us the following elegant relations. {<5j} are the roots of 

oo 

Y,Pk/(S-v k ) = l (3.10) 
fc=i 

and is given by (|3.4[) . Having found , {7^} are the roots of 

00 

X>;/0*i - "fc) = 1- (3-11) 

The last equation can be written 

A7 = 1 where A = (A kj : k,j > 1), A kj = l/(6j - v k ). 
So a formal solution is 

7 = 

Numerical solutions can be found by truncating the infinite matrix ^4 and infinite vectors 1,7 to 
N x N matrix and N- vectors, then increasing N until the desired precision is reached. 

Behaviour for large n and x independent of n. 

For p > 0, (|3.6p implies 

u n w Si/] 1 and 1^ > 0, (3.12) 

where B is given by (|3.4p . Also is given by (jA.lOp with ^ Lebesgue measure. 
Now suppose that p < 0. By (|3.9p . 

Un =71^1(1+4) (3.13) 
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where e' n — > exponentially as n — > oo and 5± has the largest magnitude among {5 n }. (The case 
where multiple 5 n exist of magnitude \8i\ requires an obvious adaptation.) 

Integral and differential equations for the eigenfunctions and resolvent. 

The right eigenfunctions satisfy UjTj = fCrj, that is 

/y 
rj((x - w)/p)dF(w) (3.14) 

For example 

Vjrj{oc) = p J rj(z)f(x — pz)dz. 

Differentiating gives the non-standard linear first order differential equation 

Vjrjiv) = sign(p)/(y)rj((x - y)/p), r,-(-oo) = 0. (3.15) 

Similarly the left eigenfunctions satisfy 
Ujlj = IjJC, that is 



ujlj(z)= / lj(y)K(y,z)dy = pf(x - pz) I lj(y)dy. (3.16) 

J Jx—pz 

So 

lj(-oo) = if p > 0, Zj(oo) = if p < 0, (3.17) 

and by differentiating, 

Uj(d/dz)[lj(z)/f(x - pz)] = p 2 l j (x - pz). (3.18) 

The resolvent satisfies 

[K{y, z, A) - K(y, z)}/\ = KK(y, z, A) = K(y, z, A)/C. (3.19) 

So 

K(-oo,z,\) = 0, 

K(y, oo,X) = 0ifp<0, 
K(y, -oo,X) = 0ifp>0 

and by differentiation the resolvent satisfies the first order partial differential equations 

(8/dy) LHSiH = S ign(p)f(y)K((x-y)/p,z,X), 
(d/dz)[LHS®n®/f(x-pz)] = p 2 K(y,x- pz,\). 

These may involve the Dirac function 5(x) since with x z = x — pz, 

(d/dy)K(y,z) = pf(x z )5(y - x z ), 

(d/dz)K(y, z) = p 2 f(x z )5(y - x z ) - p 2 I{x z < y)f(x z ). 

For special cases, it is possible to solve f)3. 14|) or ()3.16j) explicitly. 
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Example 3.1 Suppose that F(y) = ae ay on (— oo,0] where a > 0, and that Vj < 0, p < — 1, y < x. 
Taking rj(0) = 1, a solution of {3.1$ is 

rj ( y ) = e b * y where bj\vj\/a = e b > x/p , bj > 0. 

Formal expressions for the eigenfunctions. 

We now give a formal solution of (|3.15p for rj in terms of r 3 (0). (The value is arbitrary: a similar 
solution can be obtained in terms of rj(yo) for any yo.) Set 

r (y) = r j(y)i c = sign(p). 

Suppose that / and r have Taylor series expansions about 0. Denote the ith derivatives of f(y) by 
f.i(y) and set = /.i(0), n = r.j(0). Expanding 

r(y) = cf(y)r((x - y)/p) 

about 0, for i > the coefficient of y % /i\ is 

r i+1 = c^2 [ jfar.b(x/p)(-py b = c^2q ik r k = cfiro + c^2qi k r k 

a+b=i ^ ' k=0 k=l 

where 

min(i,fe) , 

m = p- k £ ( b )fi-b(-l) b x k - b /(k-b)l (3.20) 

For I, k > 1 set Qj& = qi-i yk . Set Q = (Q/ fc : l,k > 1). Set 

f = (/o, h, ■ ■ • ), R = (n, r 2 , ■ • • ), Y; = Y' = (y/1!, y 2 /2!, • • • )• (3.21) 

So R = fcr + cQR, R = ( I - cQY 1 Fcr . But r{y) - r(0) = Y'R. So we obtain the jth right 
eigenfunction in terms of its value at 0: 

r(y)/r(0) = 1 + Y^ 1 / - Q)- 1 ^ 

that is, 

r j (y)/r j (Q) = 1 + Y'((%J - Q) _1 f where ^ = v d sign(p). (3.22) 

For example for the extreme value distribution F(x) = e~ e , f = e _1 (l, 0, — 1, — 1, —7/288, —31/4, ••• )'. 

Since rj is unique only up to a constant multiplier, we may take rj(0) = 1. The solution (|3.22|) 
can now be implemented by successive approximations. For N > 1 set 

r Nj (y)/ rj (0) = 1 + Y' N (djI N - Qn)~ 1 ?n (3.23) 

where Y^v, fjy are the 1st iV elements of Y, f and Qtv is the upper left N x N elements of Q. Then 
one expects that r^jiy) —* fj{y) as N — > oo, giving the jth left eigenfunction. 

A similar treatment of (|3.18p gives an equation for the jth left eigenfunction in terms of its 
value at at an arbitrary point, taken here as x. Set 

c = p 2 \j, I = l j: e(y) = f(y)~ l . 
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By Taylor expansions, 



l(z)e(x-pz) = Y,(z i m E LH«(0)e. 6 (s)(-p) 1 



i=0 a+b= 

By (13.18p . Z satisfies (d/dz)~LHS = cl(x — pz). Taking the coefficient of z % for i > 0, 



E (*^ 1 )i.«(0)e. 6 (x)(-p) 6 = cl. i ( iC )(-p) i . (3.24) 



a+fe=i+l 

By another Taylor expansion, 

OO 



fc=0 

■VOO 



So LHS of (pOlD is £~ Wijl.jix) = Vil{x) + (WL)< where we set 



min(j',i+l) 



Wy = E e m _ a (x)(-p)^ 1 - a (-xy-7(j-a)!, 

„-n V a / 



a=0 

W = (Wij : i,i > 1), = W i0 = e. i+1 (x)(-p) i+1 , Vj = W 0j , 
Lj = l.j(x), 1/ = (Li,L 2 , ■ ■ ■), D r = diag(r { : i > 1), r = -p. (3.25) 

So (pl24"|) for i > 1 can be written UZ(sc) + W~L = cD r L so that L = (cD r - VF)~ 1 UZ(x) giving in 
the notation of (|3.2ip . 

l(z) - l(x) = Y' Z _ X L = Y' z _ x (cD r - Wy l \Jl{x). 
That is, lj(z) = l(z) is given by 

l{z)/l{x) = 1 + Y' z _ x (cD r - W)~ l \J. (3.26) 
Finally, the value of the multiplier lj(x) = l(x) is determined by ()A.3[) : 

1/ljix) = J rj {z)RHS$Fm§dz. 

An equation for the eigenvalues. 

Substituting into (|3.24p at i = 0, that is WqoI(x) + V'L = cl(x), we obtain 

V'(cD r - Wy^-U = c - Woo- (3.27) 

The roots c of this equation are just {p 2 Xj}, so this is the equation for the eigenvalues we have 
been seeking. If { 6j, j > 1} are the singular values of D~ 1 W, and if this has diagonal Jordan form 
D~ l W = R KL* where A = diag(0i,0 2 , •••), (see (fO) below), then flgjZTD can be written 

OO 

E w i( c ~ = C-00, 

i=l 

where now the weights {wi} are given by 

Wi = ViUi where v = RqD~ ~V, u = LXXJ. 
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If {cjsfj: j=i—,N+i} are the roots of its N dimensional approximation, say 

V' N {cD pN - W^Vn = c - W 00 , 

then cjyj — ► Cj = p 2 Xj as N — > oo. (This is essentially a polynomial in c of degree N + 1.) Having 
obtained an eigenvalue, one can substitute it into (|3.22p and (|3.26p to obtain the corresponding 
eigenfunctions up to constants l(x) and r(0). As noted in the appendix, either of these (but not 
both) can be arbitrarily chosen. The conditions r(— oo) = and (13.171) can be verified numerically. 

Example 3.2 Suppose that f = <j), the density of a standard normal 7V(0, 1) r.v.. Then e.j(x) = 
0(x) _1 (x) where 

H*(x) = E (x + AA(0, 1)Y = £ (ok)**'™ 7 "** 

k ^ ' 

is the modified Hermite polynomial and m,2k = (2k)\/k\2 k is the Ikth moment of M(0, 1). See 
Withers and McGavin (2006). 

An alternative is to expand RHS (|3,24p about x = 0, giving c(— p) l (l(0) + U'L) where we set 
Lj = Z.j(0), L' = (Li, L 2 , • • • )> Xk = x k /k\, U' = (xi,x 2 ,- ■ 

A ia = l. a (0)[e. b (x)(-p) b ] b=i+1 -a, A = (A ia : i,a> 1), V % = A i0 = [e. fe (x)(-p) fe ] b= , + i. 
Let 0* denote the row i— vector of zeros. For i > 1, H3.24H gives 

Vil(0) + (An, ■■■ , Ai ti+ i,0, 0, • • • )L = c{-pYiO l -\ s , xi, ■ • • )L. 

That is 

VZ(0) + = cD r XL 
where the ith row of the matrix X is (0 l_1 , xo, xi, • • • ). So 

L = S _1 VZ(0) where 5 = cD r X - A, l(y)fl{Q) = 1 + Y(,B -1 V. 

X is upper triangular, while A is lower triangular except for the 1st super-diagonal. For i = 0, 
(f3T24]) gives 

1 oo 

^Oa^a = C ZfcXfc = c/(0) + cU'L. 
a=0 fc=0 

So we obtain as an alternative equation for the eigenvalues 

A o + A i(B _1 V)i = c + cU'5 _1 V where ^ o = ~pe.i(x), A 01 = e(x). 
Unfortunately Appendix A cannot be applied with p = F since 

KG(y) = sign(p) J y G((x — w) / p)dF(w) is not of the form J K(y, z)G(z)dF(z). It would be of great 
interest, and in particular allow a unified approach to this problem, if Fredholm's theory can be 
extended to the system 

ICOr = ur, K*0*l = ul, I* Or j dp = 

for K, an q x q integral operator with kernel K(y,z) : K p x RP — > C qxq with respect a measure 
p, O a q x q operator, where * is the transpose of the complex conjugate, and v is the complex 
conjugate of v . 

For our problem, one could then apply the theory with 

P = Q = 1, V = F , K (V, z ) = sign(p)/(z < y), OG(w) = G((x - w)/p). 
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APPENDIX A 



To make the paper self-contained, we give here some theory for Fredholm integral equations 
with non-symmetric kernels. 

First consider the case where K is any k x k complex matrix. Its singular value decomposition 

is 

K = RAL* where RR* = I, LL* = I, A = diag^i, v 2 , ■■■), 
* denotes the complex conjugate transpose. Since 

KK*R = RAA*, K*KL = LA* A, 

the jth column of R is a right eigenvector of KK* with eigenvalue |Aj| 2 and the jth column of L 
is a right eigenvector of K*K with the same eigenvalue. If K is non-singular, its inverse is 

R- 1 = LA~ l R* . 

If it is singular, a pseudoinverse is given by 

K~ = LA-R*. 

However the singular value decomposition does not give a nice form for powers of K. This drawback 
is overcome by its Jordan decomposition. Consider the case where this is diagonal. Then 

K = RAL* where RL* = I, and A = diag(i/i, v 2 , ■ ■ ■ , ) (A.l) 

is composed of the eigenvalues of K. Then for any complex a, 

K a = RA a L*. 

Taking a = n and a = — 1 gives the nth power and inverse of K. 

Now let K(y,z) be a real function on 1! x where f2 is a subset of R p . Suppose that ^ is a 
cr-finite measure on $7 and that 



< \\IC\ || = J J K(y,z)K(z,y)dv(yW( 
(This L 2 condition can be changed to 



z) < oo. 



J \K(y,y)\d^y) 



< oo 



at the expense of notational complexities that need not concern us here.) The corresponding 
integral operator K, is defined by 

K4>{y) = f K(y,z)<t>(z)dfi(z), <f>{z)K = j </>(y)K(y, z)d»(y). (A.2) 

The Fredholm equations of the first kind, 

XKr(y) = r(y), \l{z)K = l(z), 
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have only a countable number of solutions, say {Xj,rj(y),lj(z), j > 1} up to arbitrary constant 
multipliers for {rj(y), j > 1}, and these satisfy 

J rjl k dfi = S jk where J gdfi = g(y)dfj,(y). (A.3) 

These are called the singular values (or eigenvalues) and right and left eigenf unctions of (K,fi) 
or IC. Also 

oo 

K(y,z)=J2r j (y)lj(z)/X j (A.4) 

3=1 

with convergence in L2(m x /x), or more strongly under other conditions: see Withers (1974, 1975, 
1978). This is the functional form of the Singular Value Decomposition for a square non-symmetric 
matrix. 

Fredholm equations of the second kind, 

r(y) - \Kr{y) = f(y), l(z) - \l(z)K = g(z), 
can be solved for A not an eigenvalue using 

(I -XIC)- 1 = I + \ICx 

where 

£\f(y) = J K(y,z,X)f(z)dfj,(z), g{z)K\ = j g(y)K(y,z,\)dfj,(y), 
and the resolvent K(y, z, X) with operator K\ is the unique solution of 

(/ - XK)K X = JC = K X (I - XIC), 

that is, 

A J K{y,u)K(u, z, X)dfj,(u) = K(y, z) — K(y, z, X) = X J K(y,u, X)K(u, z)dfi(u). 

If this can be solved analytically or numerically, then often one does not need to compute the 
eigenvalues and eigenfunctions. Alternatively, the resolvent satisfies 

oo 

K(y,z,X)=J2r j (y)lj(z)/(X j -X). (A.5) 

3=1 

The Fredholm determinant is 



D(X)=UjL 1 (l-X/X j ) = eM- J q dX J K(y,y,X)d»(y)}. 



Note that 



\\ml = Y. x i 2 - ( A - 6 ) 

3=1 
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Also since (1 - X/Xj)- 1 - 1 = A/(A j - A), 

oo 

\K( y ,z,x) = Y, r Av)h^W - VAi)" 1 - 1]. 

If only a finite number of eigenvalues are non-zero, the kernel K(y, z) is said to be degenerate. (For 
example this holds if ji puts weight only at n points.) If not, {lj} and {rj} typically both span 
L*{P) ={f-I \f\ 2 d^ < °°>- For / G L 2 (ji), 

oo oo „ „ 

f(y) = ^2 R j r j{y) = ^2 L jh(y) wh ere Rj = / //jd^x, Lj = / /r^d/i (A.7) 
with convergence in L2 (/-*)■ So 

oo oo 

K n f{y) = Y,Rjr 3 iy)/\], f{y)lC n =^Ljlj(y)/\], n > 0. (A.8) 

J=l 3=1 

So if 

| Ai| < \Xj\ (A.9) 

for j > 1 then as n — > oo, 

lC n+1 f(y)/lC n f(y) -» A7 1 , f {y)K, n+1 / f (y)K n - A^ 1 . 
This is one way to obtain Ai arbitrarily closely. Another is to use 

- 1 = sup { y gOdji : /" gr/j^t = 1} if A X > 0, (A.10) 

j" 1 = inf{ y gKhdy, : J ghdfi = 1} if Ai < 0. (A.ll) 

The maximising/minimising functions are the first eigenfunctions g = g%, h = h\. These are unique 
up to a constant multiplier if (|A.9|) holds. If Ai is known, one can use 

(Aa/C) w /(y) - flanCy), f(y)(XiJC) n -> Lx/^y), 

to approximate Rir±(y) , Lili(y) . Also since is only unique up to a multiplicative constant, 
we may choose R\ = 1 and so approximate ri(y), h(y). One may now repeat the procedure on the 
operator K\ corresponding to 

Ki(y,z) = K(y,z) - rx(y)h(z) 

to approximate X2,T2(y), assuming the next eigenvalue in magnitude, A2, has multiplicity 1. 

If say Ai has multiplicity M > 1, then 

M 

{X x KTf{y)^Y, R o r 3iyl 

3=1 

and one can adapt the method above. 

For further details see Withers and Nadarajah (In press). For further details on Fredholm 
theory for symmetric kernels, see Withers (1974, 1975, 1978). 
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